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Abstract: 



We consider the surface critical behaviour of diagonally layered Ising models on the 
square lattice where the inter-layer couplings follow some aperiodic sequence. The sur- 
face magnetisation is analytically evaluated from a simple formula derived by the diagonal 
transfer matrix method, while the surface spin-spin correlations are obtained numerically 
by a recursion method, based on the star-triangle transformation. The surface critical be- 
haviour of different aperiodic Ising models are found in accordance with the corresponding 
relevance-irrelevance criterion. For marginal sequences the critical exponents are continu- 
ously varying with the strength of aperiodicity and generally the systems follow anisotropic 
scaling at the critical point. 

PACS-numbers: 05.50.+q, 64.60.Ak, 68.35. Rh 
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I. Introduction 



The discovery of quasi-crystals[l] has stimulated intensive research to understand their 
structure and physical properties (for recent reviews see[2-6]). Theoretically a challenging 
problem is to determine the critical properties of such quasiperiodic or more generally 
aperiodic structures. After a series of numerical [7- 13] and analytical [14-21] studies on 
specific models Luck has proposed a relevance- irrelevance criterion[22]. According to this 
criterion, which is a generalisation of the Harris criterion for random magnets [23], the 
inhomogeneity is irrelevant (relevant) if the fluctuating energy in the scale of the bulk 
correlation length is smaller (greater) than the excess thermal energy. In layered systems 
the above criterion is connected to the sign of the cross-over exponent: 

0=1 + ^(0,-1) , (1) 

which is expressed in terms of the correlation length exponent v of the unperturbed system 
and the wandering exponent of the sequence a; [24]. 

Most of the studies about the critical properties of aperiodic systems are restricted 
to the quantum Ising chain with aperiodic couplings [22] [25-34], which is equivalent to the 
two-dimensional classical, layered Ising model in the extreme anisotropic limit [35]. Accord- 
ing to eq(l) for this model with v = 1 sequences with bounded (unbounded) fluctuations 
represent irrelevant (relevant) perturbations. The obtained analytical and numerical re- 
sults on different physical quantities (specific heat, surface and bulk magnetisation, local 
energy density, etc) of different aperiodic Ising quantum chains are consistent with the 
prediction of the Luck criterion. For marginal sequences non-universal critical behaviour 
was found[30-31], even if the aperiodic perturbation is of radial symmetry [34]. We note 
also on some related studies on hierarchical Ising mo dels [36]. 

In the present paper we study the surface magnetisation and the surface spin-spin 
correlations of aperiodic Ising models. In contrast to previous investigations here we con- 
sider the classical version of the model with a layered aperiodicity in the diagonal direction 
and the critical properties are studied on the (1,1) surface of a square lattice. We use two 
methods of investigations. The surface magnetisation is calculated analytically with the 
diagonal transfer matrix method[37-39], whereas both the surface magnetisation and sur- 
face correlations are numerically studied by a recursion method based on the star-triangle 
transformation. This latter method has also been used to study layered triangular systems. 
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We note that some preliminary results of our investigations have already been announced 
in a Letter [30]. 

The structure of the paper is the following. In Sections. 2 and 3 we present the methods 
of diagonal transfer matrix and that of the star-triangle recursion, respectively. Results 
on different aperiodic Ising models are given in Sec 4. A discussion is contained in Sec 5, 
while detailes of the calculation are presented in the Appendix. 



II. Surface Magnetisation by the Transfer Matrix Method 

Let us consider an Ising model on the square lattice with a diagonally layered structure, 
where the nearest neighbour couplings in the i-th layer from the surface are given by 
Ji = KiksT (see Fig la). We are interested in the magnetisation at the (1,1) surface, 
which is calculated in the transfer matrix formalism [37-39]. Denoting by < 0| and < 1| 
the ground state and the first excited state of the T diagonal transfer matrix, respectively, 
the surface magnetisation is given by the matrix-element of the surface spin-flip operator 
erf, as: 

m s =< Ojcrf jl > . (1) 

Working with free boundary conditions T is different in odd and even sites, therefore we 
consider T 2 , which is given for the inhomogeneous model as: 

AT-l 

( T2 )^7r = 2N II cosh [ K 2i-i(jH+'Pi) +K 2i (p i+1 +~p i+1 )] cosh[if 2iV -i(/UAr +~Pn)] » 
i=l 

(2) 

and depends on the configurations of the \ii = ±1 and ^ ± 1 spins, i = 1,2, . . . N (Fig. 
lb). 

To determine the eigenvectors of T 2 we make use of the fact that T 2 and the linear 
operator 

JV-l N 

H = - E A ^ 2(j *+i - E h ^ ( 3 ) 
i=i i=i 

commute 

[T 2 ,H]=0 , (4a) 
if the couplings of the inhomogeneous quantum Ising chain in eq(3) satisfy the relations: 

, C2i-2 , C 2 i+1 

h i - j ^— = h i+1 - (46) 

^2i l^2i-l 
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and 



Xi — hiS2iS2i—f 



Co,i 



2i-2 



C2i 



(4c) 



Here we used the abbreviations sinh2.ftTi = Si, cosh 2Ki = C\ and Cq = 1. Derivation 
of eqs(4a-c) is shown in the Appendix. According to eq(4a) the eigenvectors of T 2 and 
H are the same, therefore we evaluate the matrix-element in eq(l) for the inhomogeneous 
quantum Ising chain. Using a free fermionic representation of H[40] one can show[41] that: 



(5) 



and the $ s vector is determined by the equation (A + B)$ s = 0, where 

\ 

h 2 A 2 

A + B = 



(6) 



V / 

From the normalization condition £V = 1 one obtains for the surface magnetisation 



m. 



oo i / , \ 2 
hj 



1 -1/2 



i=i 7=1 vAj 



1 -1/2 



i=i j=i 



(7) 



where in the last equation we used eqs(4b) and (4c). Analysing the formula in eq(7) we 
can say that the magnetisation on the (1,1) surface of a diagonally layered Ising model 
is formally the same as that of a quantum Ising chain with inhomogeneous couplings 
Ai = sinh2Kj and in uniform transverse field hi = 1[41]. We note that for the homogeneous 
Ising model eq(7) gives Peschel's result[42]: m s = (1 — sinh -2 2K)~ 1 ^ 2 . 



III. Recursion Method 

The Ising model on the triangular lattice is invariant under the star-triangle transfor- 
mation^] (STT), which makes the exact solution of the model on this lattice relatively 
simple[44]. Also an exact renormalization group transformation for the triangular Ising 
model is based on the repeated use of the STT [45]. For a semi-infinite Ising model the 
STT has been used by Hilhorst and van Leeuwen[46] and by others[47,48] to construct an 
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iterative procedure to calculate the surface magnetisation and the surface correlations in 
the triangular Ising model. The method can be succesfully used for layered systems in 
which the couplings are the same within one layer. We note that the square lattice can 
be considered as a special case of the triangular lattice with vanishing couplings across 
the diagonals (Fig. la). In the following we shortly recapitulate the basic results of the 
recursion method. 

Let us consider a layered Ising model on a semi-infinite triangular lattice with vertical 
couplings parallel to the surface Ki, i = 1/2,3/2,... and with diagonal couplings Ki, 
i = 1, 2, . . . (Fig la). The STT maps the triangular lattice onto a hexagonal lattice which 
is in turn equivalent to a new triangular lattice. Iterating this mapping a sequence of 
triangular Ising models is generated (n = 0, 1, 2, . . .) with couplings Ki(n) and Ki(n) from 
the original model with n = 0. The surface magnetisation m s {n) and the surface spin-spin 
correlation function g s {l,n) =< aijaifi > — < o~i > 2 transform as[46]: 

m s (n) = {1 -exp [-4K 1/2 (n + 1)]} 1/2 m s (n + 1) , (8a) 

g s (l,n) =^{l-exp [-4JC 1/2 (n + 1)] } 

x [g s (l + 1, n + 1) + 2g s (l, n + 1) + g s (l - 1, n + 1)] . 

Making use of the boundary condition g s (0,n) = 1 — m 2 s {n) one obtains for the original 
model with m s = m s (0) and g s (/) = g s (l, 0)[46]: 

m s = lim [/(n)] 1 / 2 m s (n) , 

n^oo 

2n 



n=l x 



^=y:r''-[^ l jf(n)[l-ml(n)] , 



/(n) = n{l-exp[-4K 1/2 (n + l)]} . 

These relations are exact and can be used to iterate on a computer for any type of distri- 
bution of the couplings in the original layered model. In this way calculating the surface 
magnetisation one can numerically determine the T c critical point and the (3 8 critical expo- 
nent of the surface magnetisation of the model from m s (t) ~ t@ s as t = (T c — T)/T c — > 0. 
For the square lattice with Ki(0) = these results should be compared with the analytical 
expression in eq(7). 

To obtain analytical results by the recursion method one should analyse the asymptotic 
behaviour of X(n) = exp[— 4K 1 /2(n)], since according to numerical observations X(n) is 



smoothly varying with n 3> 1[48]. Inserting the asymptotic solution of X(n) into eqs(9) 
one obtains in the continuum approximation: 



r i r°° 

m s = [/(n )] 1/2 exp -- / X(n)dn 

- 1 Jna 



(10) 



and 

oo 

9s(l) = J24?2^ e M-l 2 /n)[f(n)-f(oo)] , (11) 

where no is a finite cutoff, on which the critical exponents do not depend. According to 
eq(10) the surface magnetisation is non-zero, if the integral X(n)dn is convergent, thus 
X{n) goes to zero faster than 1/n, as n tends to infinity. 

The asymptotic behaviour of X(n) has been calculated exactly at the critical point of 
the homogeneous model and for models with smoothly varying couplings at the surface [46- 
48], but there are no exact results available on X(n) outside the critical point. For the 
homogeneous, critical Ising model [46] X(n) ~ l/2n, thus m s (t = 0) = 0, f(n) ~ n -1 / 2 
and the surface correlations from eq(ll) decay as g s (l) ~ with 7711 = 1. For general 
inhomogeneous models the decay exponent follows from the asymptotic behaviour: 

lim 2nX(n) = m . (12) 

In numerical calculations it is more accurate to determine the decay exponent from eq(12), 
than to investigate the magnetisation exponent j3 s from the behaviour of the surface mag- 
netisation outside the critical point. 



IV Results on Aperiodic Models 

Although one can study general, triangular Ising models by the recursion method, 
here we restrict ourselves to the (1,1) surface of diagonally layered square models. In this 
way we reduce the space of parameters with Ki = 0, furthermore we make use of the 
analytical expression on the surface magnetisation in eq(7). 

The criticality condition for layered inhomogeneous Ising models [49] is expressed in 
terms of the variable Si = sinh 2Ki as: 

1 L 

lim - Vlog^ = (13) 

L^oo Li — 
i=l 
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Here we study two-valued sequences of the couplings and use the parametrization Si = 
Srf* , where f t takes the values or 1 according to an aperiodic sequence. The homogeneous 
model is described by r = 1. The fluctuation of the couplings in a domain of size L is 
characterised by the cumulated deviation from the average value S as [50]: 

A(S) = £(Si-S)*6L«F(^\ . (14) 

i=l ^ ^ ' 

Here S = S(r — 1) is the amplitude of the modulation, u is the wandering exponent, which is 
expressed by the leading eigenvalues of the substitutional matrix[24] to = In |A 2 |/ In Ai and 
F(x) is a fractal function of its argument with period unity. From the transformation law of 
5 under scaling one can obtain the crossover exponent 0[25] in eq(l) and the corresponding 
relevance-irrelevance criterion as described in the Introduction. The aperiodic sequences 
we consider in the following represent different types of perturbation according to this 
relevance-irrelevance criterion. 

A. Irrelevant perturbation: Thue-Morse sequence 

The binary Thue-Morse sequence[51] is generated through the substitution — > 01 
and 1 — > 10, so that one obtains after four steps: 



0110100110010110 . 



This sequence represents an irrelevant perturbation, since A 2 = and u> = — oo. 

The surface magnetisation can be obtained from eq(7) using the corresponding result 
for the Thue-Morse quantum Ising chain in Ref [27] : 



2^/2 



m, = 



y.l/2 _|_ j,— 1/2 



1 + 



1 



1 



4 V r + 1 



t + 0(t 2 



(15) 



where the critical point is at S c = r -1 / 2 and t = 1 — (S c /S) 2 . The surface magnetisation 
exponent (3 S = 1/2 takes the value for homogeneous Ising systems. Similar conclusion can 
be obtained from a study of surface critical correlations. According to numerical results 
the relation in eq(12) limn^^nX (n) =ri\\ = l is satisfied with an accuracy of 10~ 5 . Thus 
also the decay exponent takes the value for homogeneous Ising systems in two dimensions 
and the perturbation is indeed irrelevant as expected from scaling. 
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B. Relevant perturbation: Rudin- Shapiro sequence 

The Rudin-Shapiro sequence[51] is generated by the two digit substitution 00 — > 0001, 
01 — > 0010, 10 — > 1101 and 11 — > 1110, thus one obtains after three substitutions: 

000100100001110 1. 

The wandering exponent of the sequence is a; = 1/2, thus according to eq(l) this type 
of perturbation is relevant for the Ising model. The critical point from eq(13) is given 
by S c = r -1 / 2 as for the Thue-Morse model. The surface magnetisation is again obtained 
from eq(7) using the known results about the corresponding Ising quantum chain in Ref[28]. 
The surface magnetisation behaves differently for r < 1 and r > 1. For r < 1, when the 
couplings at the surface are locally stronger than the average the surface stays ordered 
at the critical point and the surface phase transition is of first order. The critical surface 
magnetisation is given by [28]: 

1 - r 

™>s,c = . r < 1 . (16) 

V 1 — r + r z 

In the other regime, r > 1, the couplings are locally weaker at the surface than in the bulk 
and the surface magnetisation behaves anomalously, it has an essential singularity at the 
critical point: 

m s ~ exp[— const (r — l) 2 t~ l ] r > 1. (17) 

According to numerical results the decay of critical surface correlations is also anomalous. 
The quantity lim n - yoo 2nX (n) — > oo, thus according to eq(12) at the critical point the 
surface correlations decay faster than any power and g 8 (l) has a stretched exponential 
dependence on /. 

C. Marginal perturbations 
Fredholm sequence 

The Fredholm sequence[51] is generated through substitution of the three letters A, 
B and C as A — > AB, B — > BC, C — > CC and we associate f t = to the letters A and C 
and fi = 1 to B. Starting with a letter A we get for the fa series after four substitutions: 

0110100010000000 . 
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This type of perturbation is localised to the surface and there is no change in the critical 
temperature, thus S c = 1. The sequence is marginal, since the corresponding wandering 
exponent u> = 0. To evaluate the formula in eq(7) for the surface magnetisation we use 
Ref [33] . For r > y/2 the surface transition is of first order, the critical surface magnetisation 
is given by: 

m s , c = ^rz| ' r>V2 . (18) 

In the other regime r < \[2 the surface transition is of second order and the corresponding 
surface magnetisation exponent is a continuous function of the parameter r: 

„ 1 lnr r- , 

Using the recursion method we have determined the decay exponent of suface correlations 
from eq(12), which is shown on Fig 2 for r < \pl. Comparing i]\\(r) with the surface 
magnetisation exponent f3 s (r) in eq(19) we can say that the surface scaling law[52] 

V\\ = 2P*/" (20) 

is satisfied for r < 

Period-doubling sequence 

The period-doubling sequence follows from the substitutional] 1 — > 10 and — > 11, 
so that starting with a 1 after four steps we have 

1011101010111011. 

The critical temperature from eq(13) is S c = r~ 2 / 3 , furthermore the sequence is marginal 
since ui = 0. The critical exponent of the surface magnetisation can be analytically deter- 
mined using eq(7) and the corresponding result for the Ising quantum chain in Ref [27]: 

= In [(l + r 2 / 3 )(l + r- 2 / 3 )] 
ls 41n2 ' 1 } 

As seen from eq(21) f3 s (r) is continuously varying with the parameter r, furthermore it is 
the same at both ends of the chain. This is a consequence of the fact that omitting the 
last digit the period-doubling sequence is symmetric. 

Next we calculate the decay exponent of critical correlations by the recursion method. 
In contrast to the surface magnetisation exponent the decay exponent is found to be 



9 



77|| = 1, independently of the inhomogeneity parameter r. On Fig. 3 we show for r = 2 
the quantity 2nX(n) as a function of the logarithm of the iterations. Its limiting value 
as n — > oo gives the decay exponent according to eq(12). The log-periodic oscillations for 
large n are a consequence of discrete scaling, which can be observed in other quantities as 
well (see eq(14)). We note that the same value of the decay exponent is found on the right 
boundary of the system. 

Comparing the surface magnetisation exponent in eq(21) and the decay exponent 
t/|| = 1 we can say that the surface scaling law in eq(20) does not satisfy. We shall come 
back to clear this point in the Discussion. 

Paper- folding sequence 

The paper-folding sequence[51] is obtained by recurrent folding of a sheet of paper, 
right over left. After unfolding one obtains a series of up- (1) and down- folds (0). The 
same sequence can be generated using the two- letter substitutions 00 — > 1000, 01 — > 1001, 
10 — > 1100 and 11 — > 1101. Starting with 11 after three substitutions the sequence is given 
by: 

1101100111001001. 

This sequence is also marginal, since u> = 0, furthermore the critical point from eq(13) is 
S c = r-V2. 

Again the surface magnetisation exponent is analytically known from eq(7) and using 
the result for the corresponding quantum Ising chain in Ref[31]. At the left surface 

= 1.(1 + r-) 
fs 21n2 V ; 



whereas at the right boundary 



_ ln(l + r) 

l: s ~ 2 In 2 ' [Z6) 



which is obtained by exchanging perturbed and unperturbed couplings, i.e. with r — > 
r _1 . Thus for the paper-folding sequence, which is not inversion symmetric, the surface 
magnetisation exponents are different at the two boundaries. 

Next we turn to calculate the decay exponent on the left boundary by the recursion 
method. Now r)\\ is found r-dependent and for all r it satisfies the relation: 

" " KW, • (24) 
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and a similar equation is true on the right surface with (3 S <-> (3 S . To illustrate the relation 
in eq(24) we show on Fig. 4 for r = 2 the quantity nX(n)(/3 s + j3 s )/f3 s , which tends to 
unity with log-periodic oscillations, in accordance with eq(12). We can say that the surface 
scaling law in eq(20) is again violated, like to the period-doubling sequence. 

V Discussion 

In this paper we have studied the surface magnetisation and the surface correlation 
function of diagonally layered Ising models on the (1,1) surface. For different aperiodic 
distribution of the diagonal couplings we have obtained exact results for the surface mag- 
netisation exponent by the diagonal transfer matrix method, whereas the decay of surface 
correlations were studied numerically by a recursion method based on the repeated use of 
the star-triangle transformation. The obtained results are in accord with the relevance- 
irrelevance criterion by Luck[22]. For the relevant Rudin-Shapiro model first-order surface 
transition and anomalous decay of critical surface correlations were observed. For marginal 
sequences (Fredholm, period-doubling and paper-folding) non-universal surface critical be- 
haviour was found, the corresponding surface magnetisation exponents are continuously 
varying with the inhomogeneity parameter r. 

The above observations remain valid, if the general triangular lattice Ising model 
with couplings K^ and Ki is concerned. Then besides the aperiodicity ratio r another 
parameter Ki/ Ki enters into the expressions. In this general case the criticality condition 
is also known analytically [49]: 



while both the surface magnetisation and the surface correlations have to be calculated 
numerically by the recursion method. Our results on the triangular lattice qualitatively 
agree with that on the diagonal square lattice, they satisfy the relevance- irrelevance cri- 
terion in eq(l) for all sequences. For marginal sequences continuously varying critical 
exponents were found, which depend on two parameters. Also the corresponding scaling 
relations are satisfied, eq(20) for the Fredholm sequence and eq(24) for the perid-doubling 
and paper- folding sequences. 

Finally, we come to the point to explain the violation of surface scaling relation in 
eq(20) for the period-doubling and paper-folding sequences. The observed scaling be- 




(25) 
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haviour in eq(24) is compatible with anisotropic scaling, when the correlation lengths par- 
allel with £|| and perpendicular to the surfaces £j_ are diverging with different exponents, 
so that £|| ~ where z is the anisotropy exponent. According to anisotropic scaling[53] 
the critical spin-spin correlation function on the left surface behave as: 

g s (l,t) = b- 2 ^g 8 (l/b z ,b 1 ^t) , (26) 

when lengths perpendicular to the surface are rescaled by a factor of b > 1. At the critical 
point t = the decay exponent is given by rj\\ — 2f3 s /vz, which corresponds to the relation 
in eq(24), if 

z = (3 s +]3 s . (27) 

For the period-doubling sequence with j3 s = f3 s , rjn = 1, as observed. We note that the 
anisotropy exponent z has been recently analytically calculated for the corresponding Ising 
quantum chains[54] in accordance with eq(27). Thus we can conclude that for marginally 
aperiodic layered Ising models where the perturbation extends over the volume of the 
system the systems become essentially anisotropic at the critical point and the anisotropy 
exponent can be expressed as the sum of the two surface magnetisation exponents. 

Acknowledgement: F.I. thanks for collaborations with L. Turban and M. Henkel 
at early stages of this work. He also thanks to I. Peschel for a useful correspondence. 
This work has been supported by the Hungarian National Research Fund under grant No 
OTKA TO12830. 
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Appendix 



To prove eq(4) we start with the representation of H in eq(3) in the (i, fx basis: 



JV-l N 



(Al) 



i=i 



i=l 



Then the matrix-elements of the commutator [T 2 , H] are given by: 



N 



cosh( J fC 2i _ 3 (^-i +7Ij_i) + K 2l _ 2 (p % -//<)) ^ _ , 



cosh(K 2i _3(^_i + + K 2i _ 2 (m + Hi)) 



i=l 



(A2) 

Here in the surface terms Kq = K-\ = K 2 n = 0. The term in the first sum in the r.h.s. 
of eq(A2) can be rewritten using the identities sinh[a( / u ± Ji)] = (fj, ± ^)/2sinh2a and 
tanh[a(// ± ~p)] = (ji ± Ji)/2 tanh2a as: 

[Hi + iHi tanh 2K 2i sinh 2K 2i - 1 cosh 2K 2i - 2 + tanh 2if 2i _ 3 sinh 2K 2i - 2 cosh 2if 2 ;-i 

+~U i+ iHi tanh 2if 2i sinh 2if 2i _i cosh 2K 2i _ 2 + Hiji^i tanh 2if 2i _ 3 sinh 2K 2i _ 2 cosh 2K 2i _i 

-(^ <-»£i)]/ 2 
so that we obtain for the commutator: 



r N-i 



I, i=l 



2 S 2 iS 2 i-i ( ft. 



2i-2 



C 2 i+i 



N-l 



C 2 i+1 



A, 



+ (Vi+lVi ~ Vi+lVi) o S-2iS2i-l (hi 3? 2 ~ h i+1 2 ' + 
7 L V °2i <^2i-l 

Then the commutator is zero if eqs(4b) and (4c) are satisfied. 



(A3) 
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Figure Captions: 



Fig. 1 (a) Diagonally layered square lattice (full line) and the corresponding triangular lattice 
with dashed vertical lines. For the square lattice the vertical couplings are zero Ki = 0. 
(b) Portion of the lattice contained in the square of the diagonal transfer matrix. 



Fig. 2 Decay exponent of critical surface spin-spin correlations for the Fredholm Ising model. 



Fig. 3 The quantity 2nX(n) as a function of the logarithm of the iterations for the period- 
doubling Ising model with r = 2. The decay exponent rj\\ = 1 in eq(12) is approached 
through log-periodic oscillations. 



Fig. 4 The quantity nX(n)((3 s + (3 S )/(3 S as a function of the logarithm of the iterations for 
the paper-folding Ising model with r = 2. The decay exponent rj\\ = (f3 s + (3 S )//3 S in 
eq(12) is approached through log-periodic oscillations. 
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